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Abstract. The light flavor antiquark distributions of the nucleon sea are calculated in the effective chiral 
quark model and compared with experimental results. The contributions of the flavor-symmetric sea- 
quark distributions and the nuclear EMC effect are taken into account to obtain the ratio of Drell-Yan 
cross sections a'^^ /2a'^'^ , which can match well with the results measured in the FermiLab E866/NuSea 
experiment. The calculated results also match the measured d(x) — u{x) from different experiments, but 
unmatch the behavior of d{x)/u{x) derived indirectly from the measurable quantity (tP°/2(tPp by the 
FermiLab E866/NuSea Collaboration at large x. We suggest to measure again d{x)/u{x) at large x from 
precision experiments with careful experimental data treatment. We also propose an alternative procedure 
for experimental data treatment. 

PACS. ll.30.Hv, 12.39.Fe, 13.60.Hb, 13.75.Cs 



1 INTRODUCTION 

The nucleon sea is an active issue in hadron physics be- 
cause of its importance in understanding both the nucleon 
structure and properties of strong interaction. In the early 
days, it was usually assumed that the sea of the proton 
was flavor symmetric between uu and dd quark-antiquark 
pairs, i.e., vP{x) = d^{x). However, this assumption was 
found to be unjustified by the observation of the Got- 
tfried sum rule (GSR) [T| violation from a number of ex- 
periments 5,6,7,8 . There have been many studies 
related to theoretical explanations of these observations 
and the flavor asymmetry of the nucleon sea 9^ . 

The Gottfried sum [T] is defined as Sq — [F2{x) — 
F2{x)]dx/x^ which, when expressed in terms of quark mo- 
mentum distributions, takes the form 



dx 



en<ifi^) + qfi^) - qn^) - €{^)]dx, (1) 



F^{x) and F^{x) 



are 



where x is the Bjorken variable, 
the proton and the neutron structure functions, is the 
charge of the quark of flavor i, and qf {qf , N=n,p) is the 
quark (antiquark) distribution of the nucleon. Decompos- 
ing the quark distribution into valance (V) and sea (Q for 
quark and A for antiquark) components, we get 



qf{x)dx 



[V^{x)+Qf{x)]dx, 



(2) 



qf{x)dx 



Afix). 



(3) 



From the flavor number conservation and the isospin sym- 
metry between the proton and the neutron, i.e., u''^{x) = 
d"(a;), mP(x) = (P{x), etc., we get 

SG = l-lj\d^x)-u^x)]dx. (4) 

If the nucleon sea is flavor symmetric, one arrives at the 
GSR 



G — o 



(5) 



The violation of the GSR was first observed by the New 
Muon Collaboration (NMC) at CERN in 1991, and 

the reported Gottfried sum is /o^[-Ff (x) — F^(x)]da;/a: = 
0.240 ±0.016. In 1994 the NMC reanalyzed their data [3], 

and got Jq[F^\x) - F^{x)]dx/x = 0.235 ± 0.026. Though 
the isospin symmetry breaking between the proton and the 
neutron at the parton level could also contribute pjOj, at 
least partially, to the GSR violation, the results are usually 
interpreted as an indication of light flavor asymmetry of 
the nucleon sea [TT], and the excess of dd pairs over uu 
pairs in the nucleon sea was measured as [3] 



[cP{x) - uP{x)]dx = 0.148 ± 0.039. 



(6) 
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The Drell-Yan process [12] can be used to measure the 
flavor distribution of the nucleon sea. The cross section of 



2 



Huiying Song et. al: Light flavor asymmetry of nucleon sea 



the Drell-Yan process at leading order is 

i 

+q^{xuQ')qf{x2,Q% (7) 

where the sum is over all quark flavors, is the charge 
of the quark of flavor i, qf- (qf^) is the quark (antiquark) 
distribution of the beam, qf {qf) is the quark (antiquark) 
distribution of the target, and xi and X2 are the Bjorken 
variables x of the partons from the beam and the target 
respectively. Two kinematic quantities commonly used to 
describe Drell-Yan process are Feynman-x, xp — xi — 
X2, and the dilepton mass M, = sa xiX2S, where 

and ^/s are the square of the invariant momentum 
transferred and the center-of-mass energy of the initial 
nucleons respectively. 

The FermiLab E772 Collaboration [13 reported an up- 
per limit on the d-u asymmetry in the range 0.04 < x < 
0.27. Later, the CERN NA51 experiment [1] measured the 
ratio of cross sections for muon pair production through 
the Drell-Yan process in pp and pD reactions at y 0, 
with 450 GeV/c incident protons. The Drell-Yan asym- 
metry is measured as Ady = (^pp - (TP")/(crPP -I- crP") = 
-0.09 ± 0.02 ± 0.025. The ratio of u over d of the nucleon 
sea derived from this measurement is 

wP 

3-|<..>=o.i8 = 0.51 ±0.04 ±0.05. (8) 

However, the acceptance of NA51 spectrometer is peaked 
near = and a; = 0.18, and consequently we are almost 
impossible to determine the a;-dependence of /vP based 
on the NA51 experiment. 

The FermiLab E866/NuSea Collaboration measured 
the ratio of cross sections of Drell-Yan muon pairs from 
800 GeV/c proton beam scattered on liquid hydrogen and 
deuterium targets and extracted d(x)/u(x) and d(x) — 
u{x) in the proton sea over a wide range of x, and got 

J^[d{x) - u{x)]dx = 0.118 ± 0.012. However, the shape 
of d'P{x)/vP{x) is beyond expectation and the ratio can 
even be less than 1 when x is large, which has received 
widespread attention because so large flavor asymmetry 
of d and u was unexpected and it is even more diffi- 
cult to explain the result that d{x) < u{x) when x is 
large. The HERMES Collaboration [S] measured charged 
hadrons from semi-inclusive deep-inelastic scattering, and 
reported d{x) ~ u{x) over the range 0.02 < cc < 0.3 and 
1 GeVVc^ < Q2 < 10 GeVVc^ Their result ofd{x)-u{x) 
is in agreement with that reported by the E866 /NuSea ex- 
periment. 

In 1977 Field and Feynman ;14 pointed out that d = u 
would not strictly hold even in the perturbative Quantum 
Chromodynamics (QCD), due to the fact that an extra 
valance up quark in the proton can lead to a suppression 
oi g ^ uu relative to g dd via Pauli blocking. Never- 
theless later calculations P^ITB] indicated that the effects 
of Pauli blocking are very small and such large asymmetry 
must have a nonperturbative origin. The role of mesons in 
DIS was first investigated by Sullivan ITT. He suggested 



that some fraction of the sea antiquark distribution of the 
nucleon may be associated with the pion cloud around the 
nucleon core. This was the original idea of the meson cloud 
model. Many authors used the pion cloud mechanism [T^ 
[TO I EO I I ^ I ^ I ^ or baryon-meson fluctuation picture 
of the nucleon to explain the light flavor asymmetry of the 
nucleon sea. Besides, the effective chiral quark model [25l 
126] is also a method to explain the nucleon sea flavor asym- 
metry [2711^ and we discuss the details in Sec. [21 There 
are many other mechanisms to explain d u, for example, 
chiral soliton model [H], instanton model [3^, statistical 
model |31j and so on. But until now, no model can explain 
the d{x)/u{x) < 1 behavior in the large x region. 

In Sec. [21 we use the effective chiral quark model US 
[26] . with the constituent quark model [32] and the light- 
cone quark-spectator-diquark model [33] results respec- 
tively as the bare constituent quark distributions inputs, 
to calculate the quark distributions of the proton. From 
the isospin symmetry between the proton and the neutron, 
we obtain the quark distributions of the neutron. In Sec. [3] 
we take the symmetric quark and antiquark sea contribu- 
tion and the Q^-evolution of quark distribution into con- 
sideration, to obtain crP'^/2crPP, d{x) — u{x) a,nd d{x) / u{x) . 
In Sec. [31 we discuss the possible nuclear EMC effect in 
the extraction of the ratio d{x)/u{x). By taking into ac- 
count such effect, the behavior of crP^2crPP, which are re- 
ally measured quantities rather than d{x) /u{x), can match 
with the experimental results better when x is large. Sec. [5] 
is devoted to some conclusions and summary. 



2 THE SEA CONTENT IN THE EFFECTIVE 
CHIRAL QUARK MODEL 

The effective chiral quark model, established by Wein- 
berg ^ , and developed by Manohar and Georgi |26j , has 
been widely recognized by the hadron physics society as 
an effective theory of QCD at the low energy scale. The 
effective chiral quark model has an apt description of its 
important degrees of freedom in terms of quarks, gluons 
and Goldstone (GS) bosons at momentum scales relating 
to hadron structure. There has been a prevailing impres- 
sion that the effective chiral quark model is successful in 
explaining the violation of GSR from a microscopic view- 
point [27112H]- Also, this model plays an important role 
in explaining the proton spin crisis |34] in Refs. |351I36) . 
A study by Ding and Xu with one of us [37] also shows 
that the strange-antistrange asymmetry within the effec- 
tive chiral quark model could explain the NuTeV anomaly. 
We adopt the effective chiral quark model to calculate the 
quark and antiquark distributions of nucleons in this pa- 
per. The principles and basic formulas are almost the same 
as those used previously, but the options and inputs are 
carefully considered and some of them are differently cho- 
sen to make the results closer to the data. 

The chiral symmetry at the high energy scale and its 
breaking at the low energy scale are the basic properties 
of QCD. Because the effect of the internal gluons is small 
in the effective chiral quark model at the low energy scale. 
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the gluonic degrees of freedom are neghgible when compar- 
ing to GS bosons and quarks. In this picture, the valence 
quarks contained in the nucleon fluctuate into quarks plus 
GS bosons, which spontaneously break chiral symmetry, 
and any low energy hadron properties should include this 
symmetry violation. The effective interaction Lagrangian 
is 



L = ^{iD^ + Vf,)Y'^ + i5A^^M7^75^ + ■ • • , (9) 



where 




(10) 



is the quark field and = 9/^ + igGf^ is the gauge- 
covariant derivative of QCD, with standing for the 
gluon field, g standing for the strong coupling constant 
and gA standing for the axial- vector coupling constant de- 
termined from the axial charge of the nucleon. and 
are the vector and the axial- vector currents which are de- 
fined by 



A.. 



1 



[11) 



where ^ = exp{in/ f), and 77 has the form 




(12) 



With the expansions for Vf^ and ^4^ in powers of 11/ f, it 
gives Vf,^0 + 0[n/ff and = id^H/f + 0{n/ff, 
where the pseudoscalar decay constant is / ~ 93 MeV. 
Thus, the effective interaction Lagrangian between GS 
bosons and quarks in the leading order becomes [17] 



(13) 



We should point out that we use the perturbative expan- 
sion in the energy rather than in the effective couple con- 
stant, which can be large. Although this model contains 
an infinite number of terms, at a given order in the en- 
ergy expansion, the low-energy theory is specified by a 
finite number of couplings. Therefore if the energy scale 
that we consider is low, the perturbative expansion is ap- 
plicable regardless of the value of the coupling constant. 
The framework that we use in this paper is based on the 
time-ordered perturbative theory in the infinite momen- 
tum frame (IMF). Because all particles are on- mass-shell 
in this frame and the factorization of the subprocess is au- 
tomatic, we neglect all possible off-mass-shell corrections. 
In this framework, we can express the quark distributions 
inside a nucleon as a convolution of a constituent quark 
distribution in a nucleon and the structure functions of a 
constituent quark. The light-front Fock decompositions of 
constituent quark wave functions have the following forms 



71' 



\U) = VZ\uo)+a,\dTT+) + ^\u7r")+aK\sK+) + ^\u7j), 

(14) 



\D) = \/Z\do) +a7r|u7r 



^\dn')+aK\sK') + ^\d^}. 

(15) 

Here, Z is the renormalization constant for the bare con- 
stituent quarks which are massive and denoted by |uo) and 
I do), and la^P are the probabilities to find GS bosons in 
the dressed constituent quark states {\U) for an up quark 
and \D) for a down quark), where a ~ Tr,K,ri. In the 
effective chiral quark model, the fluctuation of a bare con- 
stituent quark into a GS boson and a recoil bare con- 
stituent quark can be given as [55] 



qj{x) 



(16) 



In Eq. ()16p , Pja/i (?/) is the splitting function for the prob- 
ability to find a constituent quark j carrying the light- 
cone momentum fraction y together with a spectator GS 
boson a, and it has the following form 



1 / gAmV 



dk 



.y? 



(17) 



where mi^mj,ma are the masses of the i-, j-constituent 
quarks and the pseudoscalar meson a, respectively, m = 
[rrii + mj)/2 is the average mass of constituent quarks, 
and M'j^ is the square of the invariant mass of the final 
states, 



(18) 

y i^-y 

In this paper, we adopt the definition of the moment 
of the splitting function 



(19) 



with the first moment {Pja/i) — {Paj/i) = {Pa) = |aaP [55] . 
In terms of the above notation, the renormalization con- 
stant Z is given by 



Z=1-1{P^)-{Pk)-1{P 



(20) 



Now we need to specify the momentum cutoff function 
at the quark-GS boson vertex. It is conventional to use an 
exponential cutoff in IMF calculations. 



gA S^exp 



M2 



4/12 



(21) 



with — 1 following the large Nc argument [35] ■ How- 
ever, — 0.75 was adopted in the original work [55]. 
Such a form factor has the correct t and u channel sym- 
metry, and A is the cutoff parameter, which is determined 
by the experimental data of the Gottfried sum and the 
constituent quark mass inputs for the pion. This function 
satisfies the symmetry Pja/i{y) = Paj/i{^ ~ y)- 
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When probing the internal structure of GS bosons, we 
can write the process in the following form |38] 

..(x) = / ( ^) P.,,. ( ^) , (22) 

J vi V2 \yij \V2j 

where Vk/a{x) is the quark k distribution function in a 

and satisfies the normalization jj^ V}^/a{^)dx = 1. Because 
the mass of r] is so high and the coefficient is so small that 
the fluctuation of it is suppressed, the contribution of 7] is 
not considered in our calculation. From Eqs. (ITil) and (IT5|) . 
we obtain the quark distribution functions of nucleon by 
using the splitting function Eq. (1171) and the constituent 
quark distributions uq and dp, 

U{x) = Zuo{x) + Pu7v-/d (Xl do + K/7r+ ^ K+d/u <8 "O 

+ (K) (P^o„/„ (K) Uo + P^o^/a «) do), 

d(a;) = 2'do(a;) + /d7r+/M «) uo + Vd/^- » (8) do 

+ ^Pd-RO/d «> do + Vd/Ko «) Pxos/d « do 

+ ^Vd/„o «) (P^o„/„ «) Uo + P^od/d «) do)- (23) 

Here, we define the notations for the convolution integral 
as 

dy 



P 



® * = / — ^ja/» (2/) q^ (-) , (24) 



and 



r ^ {If''- fe)-- (I) 

In the same way, we can derive the light-flavor antiquark 
distributions, 

u{x) = Vf,/^- (g) P^-u/d ® do 

+ ^Va/^o «) (P^o„/„ (8) Uo + PrrOd/d do), 

d(a;) = Vj/7r+ «) P^+<i/« «) "0 

+ ^Kj/ttO (PttOu/u <8) Uo + P^od/d <E) do). 



(26) 



where 



and 



Kt/7r+ = V'd/7r+ = ^d/TT ^fi/TT- 

= 2K/7I-0 = 2t4/7rO = 2Vd/^o = 2T/j/^o 

= (a;) , 



Vu/K+ = Vd/K" ■ 



(27) 



From above equations, we can reexamine the valence quark 
distributions Uv{x) — u{x) — u{x) and dy{x) = d{x) — d(x), 
which satisfy the correct normalization with the renormal- 
ization constant Z. We should point out that in the chiral 
quark model, antiquarks are produced in the process of the 
splitting of Goldstone bosons unless higher order effects 
are considered. Since Goldstone bosons are spin-0 parti- 
cles without polarization, antiquarks are unpolarized, and 
this feature is compatible with the available experimental 
data [40] . 

In this paper, we choose m„ — — 330 MeV, mg = 
480 MeV, m^± — m^o — 140 MeV and m/j-+ — mj^o = 
495 MeV. Employing the quark distributions of the ef- 
fective chiral quark model, we get the Gottfried sum de- 
termined by the difference between the proton and the 
neutron structure functions, 



Sg= f'-[F^{x)-F-ix)] 
Jo X 



da;[u(a;) -I- u{x) — d{x) — d{x)] 



= liZ^l{P,) + {PK) + ^{P,)) 

-^(1-2 (P.)). 



(28) 



From the above equation and the experimental value of 
Gottfried sum [3 , we can find that the appropriate value 
for is 1500 MeV. At the same time, (P^) = 0.149, 
{Pk) = 0.085, (P^) = 0.063 and Z = 0.682. But for K and 
f] mesons, the terms (Pk) and (P^) in the Gottfried sum 
cancel out those terms in Z = 1 — ^{Ptt) — {Pk) — ^{P^)^ 
so the value of A can not be determined from Eq. 
or experimental data. It is natural to assume that the 
cutoffs are same for tt, K and rj mesons in the effective 
chiral quark model, A^, = Ak = 1500 MeV 38,41, which 
is different from the traditional meson cloud model. 

In this paper, the parton distributions for mesons are 
taken from the parametrization from GRS98 given by Gluck- 
Reya-Stratmann [42j because the results from parametriza- 
tion can be closer to reality than those from models. 



V^(x) = 0.942X 



-0.501/ 



H-0.632^/i)(l 



\0.367 



K/K+(a;) = Fd/A'o(x) = 0.541(1 - xf-"V^{x). (29) 

We also need inputs of constituent-quark distributions 
Uo and do. But there is no proper parametrization of them 
because they cannot be directly measured in the experi- 
ment. Therefore, we have to choose some models as in- 
puts. In this paper the constituent quark model distri- 
butions and the light-cone quark-spectator-diquark 
model distributions 33 are adopted as two different kinds 
of inputs of constituent quark distributions. The constituent 
quark model distributions have the following forms 



uo(a;) 



-x^-^iX^x) 



C1+C2 + 1 



B[ci + l,ci +C2 + 2]' 

'^o(^) = Br ^/o ^., x^\\-xf^-^\ (30) 
B[C2 + 1, 2ci + 2J 
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where B[i,j] is the Euler beta function, and ci — 0.65 
and C2 = 0.35 adopted from Ref. [32ll43] can satisfy the 
number sum rules 



1 .1 
uo{x)dx = 2, / do{x)dx = 1, (31) 

^0 



and the momentum sum rule 



1 .1 
xuQ{x)dx + / xdQ{x)dx = l. (32) 



It is pointed out that there are other different values for 
ci and C2 suggested by Ref. [35|. The light-cone quark- 
spectator-diquark model distributions [33] are 

uo{x) = 7;as{x) + ^av{x), 

2 D 

do{x) = ^av{x). (33) 



where aoix) cx J[d'^k±]\ipDix,'k±)\'^ {D — S or V) is 

normalized such that dxauix) = 3 and denotes the 
amplitude for the quark q being scattered while the spec- 
tator is in the diquark state D. We adopt the Brodsky- 
Huang-Lepage prescription |45) for the light-cone momen- 
tum space wave function of the quark-spectator-diquark 



LpD{x,'k±) = ^Dexp{- 



1 



l-x 



-]}• 

(34) 

where k_L is the internal quark transversal momentum, 
rUq and m^i are the masses of the quark q and specta- 
tor D, and is the harmonic oscillator scale param- 
eter. In this paper we simply adopt niq — 330 MeV, 
/3d = 330 MeV, ms = 600 MeV and my = 800 MeV 
as often adopted in literature. In the light-cone quark- 
spectator-diquark model, the number sum rule Eq. (j3ip 
is still satisfied but the momentum sum rule Eq. (1321) is 
violated. It should be noticed that in this paper we adopt 
the quark-spectator-diquark model, i.e., when any quark 
is probed, the other part of the target is served as a spec- 
tator with quantum numbers of a diquark. Thus, some 
gluon effects may exist inside the spectators. This means 
that partial momentum can be carried by gluons at the 
initial point within the quark-spectator-diquark model. 
Hence there is no need to require the momentum sum 
rule Eq. ([5^ by quarks as in the constituent quark model, 
where the nucleon momentum is distributed among con- 
stituent quarks at the initial point. 



3 ADDITIONAL SYMMETRIC SEA 
CONTRIBUTIONS 

As is well-known, the quark distributions measured by ex- 
periments at certain include not only non-perturbative 
intrinsic sea but also perturbative extrinsic sea [53]. Al- 
though the antiquark content in the nucleon mainly comes 



from the intrinsic sea, the extrinsic sea should also be con- 
sidered from a strict sense when we want to investigate the 
distributions of quarks and antiquarks. Therefore, we take 
into account additional contribution from the symmetric 
sea before we use our quark distributions to compare with 
experimental data. 

We should point out that the quark distribution func- 
tions we obtain in the front can only be proper at a cer- 
tain Qq, because the values of parameters in the model 
do not evolve according to . In the experiment of the 
E866/NuSea Collaboration, varies from about 21 GeV^/c^ 
to more than 160 GeV^/c^, so the change of should 
be considered carefully. In this paper, we choose Qo = 
0.7 GeV/c. For simplicity, we also assume that the quark 
distribution functions we get have the same evolution be- 
havior with that of parametrization we adopt. 



q (x,g^) = q {x,Q, 



2N 9 




pare 



para 



{x,> 



(35) 



where q = d, u, d. The gP"** stands for the quark dis- 
tribution of flavor q from parametrization. 

It is found that most parametrizations of parton dis- 
tributions after the experiment of E866/NuSea Collabo- 
ration had been considerably affected by it, so we adopt 
an earlier parametrization, namely CTEQ4 parametriza- 
tion |46j . In this paper we will use u as the criterion. Other 
methods to take the contribution of symmetric sea may 
also work, and our method is a feasible choice. As we can 
see, the u{x) we derived from the model can be larger 
than that of CTEQ4 parameterization when x is large, so 
we should not add symmetric sea contribution to u and 
d any more when they are larger than the result from 
parameterization. Specifically, if ■it™°'^°'(a;) < uP'^'''^{x), we 
estimate the symmetric sea contribution du{x) = Sd{x) 



Su{x) = 5d{x) = uP"^(x) 



odd 



(x), otherwise we set 



the symmetric sea to be zero. Intuitively, the flavor sym- 
metric sea perturbative extrinsic can be thought as arising 
from the splitting of gluons to quark- antiquark pairs, thus 
they are flavor symmetric. Similar consideration has also 
been adopted to confront the calculated strange and an- 
tistrange distributions with experimental observations in 
the calculation of the strange- antistrange asymmetry in 
the chiral quark model [37) . 

We compare the parton distributions derived above 
with the CTEQ4 parametrization and find that u and d 
quark distributions can match well while that u and d 
quark distributions really vary. Thus, we will adopt two 
methods to compare results with experimental data in the 
following: (1) we use parton distributions of both quarks 
and antiquarks from the model in the calculation, and this 
is denoted as "Model" ; (2) we use parton distributions of 
u and d from the model while the distributions of u and 
d are parametrization results from CTEQ4 directly, and 
this is denoted as "Q.Para.". 

The forms of cross sections we take are 

4 _ 1 - 

cx -u{xi)u{x2) + -d{xi)d{x2) 

4_ 1 - 

-f -u{xi)u{x2) + -d{xi)d{x2), (36) 
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Fig. 1. Distributions for xu{x), with = 0.49 GeV^/c^ for the left part and = 54 GeV^/c^ for the right part. The dotted 
curve is the result from CTEQ4 parametrization. The thin solid and dashed curves are the model calculation results with the 
constituent quark model (CQM) and the light-cone quark-spectator-diquark model (LQM) as inputs. The thick corresponding 
curves are the corresponding results with chiral quark model results plus the symmetric sea contributions. 




Fig. 2. Distributions for xd{x), with = 0.49 GeV^/c^ for the left part and = 54 GeVVc^ for the right part. The dotted 
curve is the result from CTEQ4 parametrization. The thin solid and dashed curves are the model calculation results with the 
constituent quark model (CQM) and the light-cone quark-spectator-diquark model (LQM) as inputs. The thick corresponding 
curves are the corresponding results with chiral quark model results plus the symmetric sea contributions. 



and 

4 - 1 _ 

ctP" cx -u{xi)d{x2) + -d{xi)u{x2) 

4_ 1 - 

+ -u{xi)d{x2) + -d{xi)u{x2). (37) 

The influence of heavier quarks is not included as their 
contributions can be reasonably neglected. With all these 
taken into consideration, we can get the results shown in 
Fig.[3]and Fig. 01 We find the results are quite good except 
when X is large, regardless of the constituent quark dis- 
tribution inputs and the methods we adopt. Accordingly, 
the approach taken by us should be reasonable. 



4 DISCUSSION ON NUCLEAR EMC 
EFFECT 

In 1982, the European Muon Collaboration (EMC) at 
CERN |47','48; found that the structure-function ratio of 
bound nucleon to free nucleon, in the form of F^{x, Q^)/ 
F^{x,Q'^), is not consistent with the expectation by as- 
suming that a nuclei is composed by almost free nucle- 
ons with Fermi motion correction taken into account, and 
such phenomenon was confirmed by E139 collaboration 
at SLAC (49j. This discovery, which is called the nuclear 
EMC effect, has received extensive attention by the nu- 
clear and hadronic physics society. There are many mod- 
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Fig. 3. The Drell-Yan cross section ratio a^^ /2a^^ versus X2. 
The data are taken from E866/NuSea [7] experiment. The 
following curves are results from model with the light-cone 
quark-spectator-diquark model (LQM) as inputs plus symmet- 
ric sea contributions. The thin dash-dotted (dotted) curve is 
the result that the quark distribution from the chiral quark 
model (CTEQ4 parametrization) without EMC effect. The 
thick dashed (solid) curve is the corresponding result with 
EMC effect for the parameter ^ = 1.1. The thin dashed (solid) 
curve is the corresponding result with EMC effect for the pa- 
rameter ^ — 1.5. 
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els describing EMC effect now ISOlEIlEllSllMlESlEaEZl 
[58] , and a good review can be found in Ref. [59] . All these 
models can qualitatively describe the data in the medi- 
ate X region. The inclusive deep inelastic scattering (DIS) 
data are expressed as F^{x, Q^), which can be 

written in the naive parton model as 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 



Fig. 4. The Drell-Yan cross section ratio a'^^ /2a'^'^ versus 
X2- The data are taken from E866/NuSea |7j experiment. The 
following curves are results from model with the constituent 
quark model (CQM) as inputs plus symmetric sea contribu- 
tions. The thin dash-dotted (dotted) curve is the result that 
the quark distribution from the chiral quark model (CTEQ4 
parametrization) without EMC effect. The thick dashed (solid) 
curve is the corresponding result with EMC effect for the pa- 
rameter ^ = 1.1. The thin dashed (solid) curve is the corre- 
sponding result with EMC effect for the parameter ^ = 1.5. 



F^{x,Q^) _ S^e| [g,(x,g^A)+g,(x,g^A)] 



i^D(a:,g2) 



r,ef Mx,Q^) + q,{x,Q^)] 



(38) 



where denotes the charge of the partons with flavor 
i, and q{x,Q^) is the parton distribution function of the 
nucleon. 

In the analysis of the E866/NuSea Collaboration, the 
nuclear effects in deuterium were assumed to be negligi- 
ble. As discussed above, the calculated results oia^^ /2a^P 
cannot describe the experimental data well when x is large, 
and this is just the region where the nuclear EMC effect 
may begin to work. So we will take the EMC effect into ac- 
count to check how the results can be changed. In this pa- 
per, we choose the so-called Q^-rescaling model [5^15711551 
[55] . In this model, the quark of a bound nucleon in the nu- 
clear medium is considered to have different confinement 
size compared with that of the quark in the free nucleon, 
and consequently q^{x, Q^) is related with q^{x, Q^) (the 
parton distribution in the free nucleon) by the relation 



q^ (x,g2)=<z^ {x,aQ') Q') 



(39) 



At Q Ri 12.5 GeV/c, we adopt ^ = 1.1 given in the orig- 
inal work • We show the ratio of (x) in deuterium 
and that in a free proton plus a free neutron (denoted as 
p-|-n) in FIGlSJ We assume furthermore that the nuclear 
EMC effect only takes effect when x > 0.22. To show 
the dependence of the rescaling factor ^, we also adopt a 
larger value ^ = 1.5 as a comparison. Then we can see 
that the behavior of cross section ratio at large x is visi- 
bly improved. In addition, we find that the cross section 
rjitio at {X2) — 0.315 is smaller than 1. We also display 
d{x) — u(x) and d{x)/u{x) in FiglHland Fig]?] respectively. 
As we can see, the behavior of d(x) — u{x) we derived can 
match well with the experimental data and parametriza- 
tion from CTEQ4, however, d{x)/u{x) we get can match 
with experiment at small x but is not compatible with the 
experiment at large x, and especially we could not get the 
ratio smaller than 1 in the region we consider. It is wor- 
thy to remind that all known models have the result that 
d{x)/u{x) is larger than 1 over all x range. 

It is found that we can adopt an alternative proce- 
dure to deal with u quark and d quark to compare the 
result with the experiment. We assume u{x) = uP'^'^^{x) 
and d{x) = dy{x) + d{x), where d^ = d^'^'"'^ - Jp'*'''' is the 
valance d quark distribution, and d is the result we get 
from the chiral quark model with symmetric sea quark 
distribution added. The result is denoted as "Dv-Para." 
and displayed in Fig. |SJ We find that this result can also 
match well with the experimental data. 

So we suggest a procedure for experimental data treat- 
ment to derive d{x)/u{x) and d{x) — u{x) from the cross 
section ratio a'^^ /2a'^^: first assume u{x) ~ vP^^^{x)^ u{x) - 
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X 

Fig. 5. The ratio of bound deuterium F^{x) compared to free 
F2^'^{x). The sohd (dashed) curve is the resuh with ^ — 1.1 
(C = 1.5). 



^paraj-^) and d{x) = dP'^'^^ix) - dP^'^^ix) + d(x),_then the 
experimental data is used to fit the behavior of d{x). This 
procedure might give results more compatible vifith model 
predictions than the method adopted in the E866/NuSea 
analysis. 



5 RESULT AND CONCLUSION 

In this work, we calculate the light flavor quark and anti- 
quark distributions within the effective chiral quark model 
by using the constituent quark model and the light-cone 
quark-spectator-diquark model as inputs respectively, and 
revise the results by taking into consideration the sym- 
metric nucleon sea contributions, the Q^-evolution and 
the nuclear EMC effect. The distributions of (tP°/2ctPP 
and d{x) — ■u{x) match with the experimental data, and 
d{x)/u{x) is compatible with the experiment at small x, 
while the behavior of d{x)/u{x) in large x region is differ- 
ent from the experimental result. However, the result di- 
rectly measured in E866 experiment was only (tP^/2(tPp, 
whereas d{x)/u{x) was derived indirectly from crP'^/2crPP 
with several assumptions. Although it is entirely possible 
that the analysis of E866/NuSea Collaboration was cor- 
rect, it is also possible that the analysis of them was based 
on some excessive assumptions. For example, the assump- 
tion of u{x) + d{x) fixed as parametrization may have a 
strong influence on the ratio d{x)/u{x) derived from data 
of cross sections. In addition, it is worth noting that the 
nuclear EMC effect should be considered carefully when x 
is large. Therefore, without ruling out these possibilities 
and considering all we discussed before carefully, the re- 
sults of d{x)/u{x) extracted from other quantities might 
be not so reliable as the results of crP^/2crPP. We also 




-0 2 I I I I I I I I I I I I I I I 

0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 

X 

Fig. 6. The distribution for d{x) — u{x). The data are from 
HERMES m and E866/NuSea [7] experiments. The sohd and 
dashed curves are the model calculation results in the effective 
chiral quark model with the constituent quark model (CQM) 
and the light-cone quark-spectator-diquark model (LQM) as 
bare quark distribution inputs respectively. The dotted curve 
is the result from CTEQ4 parametrization. All the values are 
scaled to fixed = 54 GeVVc^ 



2.4 




2 I I I I I I I I I I I I I I Li I 

0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 



X 

Fig. 7. The distribution for d{x)/u{x). The data are from 
NA51 g| and E866/NuSea [7] experiments. The solid and 
dashed curves are the model calculation results in the effective 
chiral quark model with the constituent quark model (CQM) 
and the light-cone quark-spectator-diquark model (LQM) as 
bare quark distribution inputs plus symmetric sea contribu- 
tions. The dotted curve is the result from CTEQ4 parametriza- 
tion. All the values are scaled to fixed = 54 GeVVc^. 
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Fig. 8. The DreU-Yan cross section ratio a'^^ /2a'^'^ versus X2. 
The data are taken from E866/NuSea [7] experiment. u{x), 
u{x), dvix) are obtained by adding CTEQ4 parametrization 
and symmetric sea contributions. The thin dash-dotted and 
dotted curves are the results with d{x) derived from the chiral 
quark model with constituent quark model (CQM) and light- 
cone quark-spectator-diquark model (LQM) as bare quark dis- 
tribution inputs without EMC effect. The thick solid and 
dashed curve are the corresponding results with EMC effect 
for the parameter ^ = l.L The thin solid and dashed curves 
are the corresponding results with EMC effect for the param- 
eter ^ — L5. 



suggest an alternative procedure to derive d{x)/u{x) and 
d{x) — u{x) from the experimental data of crP^/2o'PP. Thus 
it is important that more precision experiments should be 
carried out to enable more direct and accurate determi- 
nation of sea quark and antiquark distributions. 
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